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Anyons from non-solvable finite groups 
are sufficient for universal quantum computation 

Carlos MochorQ 
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We present a constructive proof that anyonic magnetic charges with fluxes in a non-solvable finite 
group can perform universal quantum computations. The gates are built out of the elementary 
operations of braiding, fusion, and vacuum pair creation, supplemented by a reservoir of ancillas 
of known flux. Procedures for building the ancilla reservoir and for correcting leakage are also 
described. Finally, a universal qudit gate-set, which is ideally suited for anyons, is presented. The 
gate-set consists of classical computation supplemented by measurements of the X operator. 

PACS numbers: 03.67.Lx, 05.30.Pr 



I. INTRODUCTION 

The discovery of the potential speedups offered by 
quantum computers launched an effort to find physical 
systems out of which these computers could be built. 
Researchers soon found that these systems are in short 
supply, as it is extremely difficult to isolate a quantum 
system from the environment, while maintaining enough 
control to perform operations on the encoded data. The 
advent of quantum error correction and fault-tolerant 
processing has drastically increased the tolerable error 
rates; nonetheless, physical systems with low enough er- 
ror rates are still hard to come by. 

One way to protect a Hilbert space from the environ- 
ment is to encode the quantum data in non-local observ- 
ables. These observables, which are constructed from 
topological invariants, cannot be measured or changed by 
any local operator. Therefore, because the environment 
generally acts locally, the physics of the system provides 
a form of fault tolerance. 

In particular, consider the spectrum of electrically 
and magnetically charged particles that are obtained by 
breaking a gauge group to a finite subgroup. The fi- 
nite group gauge theory is a particularly good system for 
quantum computing because it involves no gauge fields, 
and hence no long-range interactions except for those ob- 
tained by braiding. Furthermore, the Hamiltonian of 
the system respects the unbroken symmetry; therefore 
Schur's lemma forbids the types of coupling to uncharged 
objects that can produce decoherence. Of course, the 
data could still decay by the exchange of a charged par- 
ticle between two anyons, but this is a quantum tunneling 
event which is exponentially suppressed by the distance 
between the particles. 

When the gauge theory is restricted to two spatial 
dimensions, the particles acquire topological long-range 
interactions, which can be be used to perform compu- 
tations. The interactions occur when the particles are 
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exchanged or braided, and depend only on the topologi- 
cal class of the path involved. Because of these interac- 
tions, the charged particles have quantum statistics that 
are more exotic than the standard fermions and bosons, 
and are known as anyons. The non-standard statistics, 
though, only arise when clockwise rotations can be distin- 
guished from counterclockwise rotations, which is why we 
impose the requirement of two spatial dimensions. While 
this two-dimensional model of the world seems somehow 
unphysical, there exist condensed-matter systems with 
quasiparticles that behave like anyons. 

The original proposal for an anyon based quantum 
computer was made by Kitaev Q, 0- The first con- 
crete description was done by Ogburn and Preskill in 
Refs. 0,3 for anyons in the group A$, the even permu- 
tations of five elements. In our paper, we will generalize 
the work of Ogburn and Preskill to any non-solvable fi- 
nite group, which includes A$ as the smallest case. 

The paper is organized as follows: We begin by intro- 
ducing some notation and reviewing the properties of the 
anyon model that will be used throughout this paper. 
Section II I II presents the universal gate-set that will be 
employed to prove anyons can perform quantum compu- 
tations. Sections IIVI through IVII contain the meat of the 
paper, and discuss a concrete anyonic implementation of 
all the necessary gates. For pedagogical reasons, we first 
cover the easier subcase of simple perfect groups in Sec- 
tion llVI and then discuss the required generalizations for 
any non-solvable group in Section IVII In Section we 
discuss how to make these computations fault-tolerant 
by performing leakage correction. Finally, we discuss the 
conclusions and unsolved questions. There are also two 
appendices which include the mathematical proofs, and 
a technique for creating anyon ancillas. 



II. REVIEW 

In this section we will review some of the braiding and 
fusion properties of our anyons. Our review will be rather 
abridged, but more details can be found in the excel- 
lent review of discrete gauge theories @ (and the original 
work JjJ). The paper by Ogburn and Preskill 0,0] also 
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FIG. 1: Exchanging two anyons. 



contains a good review with emphasis on the applications 
to quantum computing. 

This section also establishes our notation for qudits, 
and reviews the phase estimation circuit, a highly useful 
trick that will be used often. 



A. Magnetic charges 

The main players throughout this paper will be the 
magnetic charges, also known as fluxes. For a field the- 
ory with an unbroken finite group G, there is one mag- 
netic charge for each element g 6 G. Quantum mechan- 
ically, we can have superpositions of these states, giving 
a one-particle Hilbert space spanned by \g) for all g G G 
(though strictly speaking, superpositions of charges in 
different conjugacy classes are meaningless, as will be 
explained in the next subsection). 

Specifying the exchange properties of the charges in- 
volves making a choice of gauge. The easiest choice, 
which will be used in this paper, is to keep all anyons 
ordered on a horizontal line. The exchange of particles, 
which can be clockwise or counterclockwise, is only al- 
lowed between adjacent pairs. In either case, the particle 
that passes below remains unchanged, while the particle 
that passes above gets conjugated. When the exchange is 
in the counterclockwise direction, the upper anyon gets 
conjugated by the flux of the lower one, whereas in the 
clockwise direction it gets conjugated by the inverse of 
the lower flux. This is depicted in Figure ^ 

One way to visualize these exchanges is to associate 
with each anyon a ray that is vertical in the plane, start- 
ing at the particle and proceeding upwards. Anyons are 
allowed to move freely through the plane, but every time 
an anyon crosses the ray of another particle, it gets conju- 
gated by the flux of the owner of the ray (or by the inverse 
flux if crossing from left to right). Note that when a par- 
ticle passes a group of anyons, it gets conjugated by the 
total flux of the anyons, which is given as the product 
from left to right of the individual fluxes. 

Clearly, moving single anyons around can produce 
strange correlations throughout the system. However, 
moving a pair of anyons with a total flux that is trivial 
will not change the state of the system if the pair always 
passes below. This is why we will always be dealing with 



FIG. 2: Conjugating a pair of anyons. 
states of the form 

E^l^lT 1 }, (1) 
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which correspond to a pair of anyons with trivial total 
flux. When dealing only with pairs of trivial total flux, we 
can swap any two pairs, or bring any two pairs together 
without affecting the state of the rest of the system. 

We do want to allow controlled interaction between 
pairs, though, and this is accomplished by a pass-through 
operation. The idea is to have one pair circle one anyon 
from the other pair. This will conjugate the fluxes of 
the pair that circles, but leave the other pair invariant. 
This operation is depicted using elementary exchanges in 
Figure |21 

The net result of the pictured operation is 

\h) <8> Ih- 1 ) ® \g) ® Ig- 1 ) 

— > \h) ® Ih' 1 ) g> \hgh~ 1 ) <g> Ihg^h' 1 ), (2) 

which is a conjugation of the second pair by h. Conjuga- 
tion by h^ 1 could be achieved by using counterclockwise 
exchanges in the picture. 

For notational convenience, in this paper we will gen- 
erally only mention the flux of the left element of a pair. 
The above transformation will be written as 

\h)®\g)^\h)®\hgh- 1 ), (3) 

leaving the compensating fluxes implicit. While we will 
exclusively deal in this paper with flux pairs with trivial 
flux, we will only explicitly refer to the second anyon 
when necessary to describe the operations. 
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B. Electric charges and vacuum pairs 

We now wish to focus on the operations of creating 
pairs from the vacuum and fusing pairs back into the vac- 
uum. However, we must first briefly discuss the complete 
spectrum of particles, and that involves electric charges. 

An electric charge is a particle with no flux that trans- 
forms as some non-trivial irreducible representation of 
the group G. A useful analogy is to think of the rep- 
resentation of G as the total spin of the particle. The 
internal state of the particle is then equivalent to the 
direction in which the spin is pointing. 

The electric charge states can be labeled as \R,V), 
where R is a representation of G and V is a vector that 
transforms in the representation R. The electric charges 
do not interact with each other, but when one of them 
circles a magnetic flux g, its state changes as 

\R,V) — \R,U R (g)V), (4) 

where Un(g) is the matrix corresponding to g in the rep- 
resentation R. This is known as the Aharonov-Bohm 
effect. 

While we can transform the state of an electric charge 
within the subspace of a representation, there are no op- 
erations (other than fusion, which destroys the particle) 
that can change the representation of a particle. Fur- 
thermore, the phase between states of different represen- 
tations cannot be measured. We can therefore effectively 
describe the electric charges as having decohered into the 
different representations. In particle physics we would 
say that the different representations correspond to dif- 
ferent superselection sectors. 

The same thing happens to the magnetic charges. Dif- 
ferent conjugacy classes live in different superselection 
sectors, so we can imagine that there is an automatic 
decoherence into different conjugacy classes. Superposi- 
tions of fluxes in different conjugacy classes are therefore 
meaningless. 

The spectrum also contains particles with both elec- 
tric and magnetic charge, which are called dyons. The 
only special feature is that the electric charge is a repre- 
sentation only of the subgroup of G that commutes with 
the flux. The aforementioned magnetic charges are just 
dyons with a trivial representation. The dyons also have 
superselection sectors that correspond to different conju- 
gacy classes and representations. 

The purpose of discussing the full spectrum, and the 
idea of superselection sectors, is to find out what kind 
of states we get when we create a pair of particles from 
the vacuum. The first thing to note is that each of the 
particles will instantly decohere into a specific conjugacy 
class and representation. Furthermore, because a pair 
created from the vacuum must have trivial total charge 
and flux, the conjugacy classes must be inverses, and the 
representations must be conjugate representations. 

Consider the case that the pair decoheres into plain 
magnetic charges, with the first one contained in the con- 
jugacy class C. Because the combined state still has vac- 



uum quantum numbers, the state must transform triv- 
ially when another flux is dragged around it. That is, it 
must be invariant under conjugation. There is only one 
such state: 

iVac^H-^^Tl^kT 1 ). (5) 

VIM g ec 

The vacuum states for the other superselection sectors 
are also unique and have similar expressions. When a 
pair of anyons is created from the vacuum, it will start 
off in one of these states. 

Another useful operation is to fuse two anyons to- 
gether. Note that we are not talking about two anyon 
pairs, but rather two anyons, sometimes from the same 
pair, and sometimes from different pairs. The operation 
of fusion will turn the two particles into one, which must 
carry the total flux and charge of the two. It is also 
possible that the two anyons will have vacuum quantum 
numbers, and will fuse back into the vacuum. In this 
case, no particle will be left behind and their total mass 
will be transformed into some other medium, such as ra- 
diation. If is the combined state of the two anyons, 
and the first anyon is in the conjugacy class C, then the 
probability that the two will fuse into the vacuum is given 
by the standard rules of quantum mechanics: 

Pvacuum — 

|(Vac(C)|vI/)| 2 . (6) 

After fusing two particles of different pairs, the fused 
particle may carry some flux. However, since the total 
flux of the original four particles was trivial, the total 
flux of all the remaining particles (including the product 
of the fusion) will be trivial as well. Therefore, it is pos- 
sible to safely move the group of particles away from the 
bulk of the computation without disturbing our quantum 
state. 



C. Qudits 

Throughout this paper it will be useful to perform 
computations with qudits rather than the usual qubits. 
We define our computational basis as the states \i) for 
< i < d, where we will assume that d is prime. The 
unitary Z and X gates can be defined as follows 

Z\i) = w*!*), (7) 
X\i) = \i + l), (8) 

where u> is a fixed non-trivial d th root of unity, and sums 
are understood to be modulo d. The operators satisfy 
the commutation relation 

ZX = XZuj. (9) 

As usual, the eigenstates of Z correspond to the com- 
putational basis. We can also introduce the eigenstates 
of X: 
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which have the following transformations under the ac- 
tion of our unitary gates: 



Z\i) = 
X\i) = 

D. Phase measurement 



(11) 
(12) 



A very useful trick, used many times throughout this 
paper, is Kitaev's phase estimation technique In fact, 
we will only employ a special case of the technique which 
we describe below. 

Assume that we are working in a system with qudits, 
and we have an operator U with eigenvalues that are d th 
roots of unity. We shall prove that being able to apply a 
controlled-?/, and measure in the X basis, is equivalent 
to being able to measure the operator U, 

Consider applying the circuit below to an eigenstate 
of U with eigenvalue uP: 
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u- 1 









* + 3) 



I*, 



where the controlled-?/ -1 can be applied as d — 1 
controlled-C/'s. The circuit works as described because 
the controlled- C/ _1 leaves the bottom state invariant, but 
applies a Z~ J to the upper state. On a general state 
\4>) = J2j c j\^j) expanded in terms of eigenvectors of U, 
the circuit produces the transformation 



(13) 



Clearly, a subsequent measurement of the first qudit in 
the X basis is equivalent to a non-destructive measure- 
ment of the original state in the U basis. We will use 
this technique in the next section to measure the opera- 
tors X a Z b . 

In a later section we will employ the equivalent circuit 
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x- 1 









run in both the forwards and backwards direction, to 
change between the \i) states and the readily available |0) 
state which can be naturally produced from, and fused 
into the vacuum. 



III. A UNIVERSAL GATE-SET FOR ANYONS 

A lot of the work in proving universality can be simpli- 
fied by choosing a proper gate-set. For this paper, we will 
employ a generalization of the gate-set used by Ogburn 
and Preskill 0,3- The gate-set, which involves measure- 
ments as well as unitary gates, can be applied to qudits 
when d is prime, which is the only case considered in this 
paper. 

The universal gate-set is: 

1. Measure non-destructively Z 

2. Measure non-destructively X 

3. Apply Toffoli operators (to any set of three qudits) 
where the qudit Toffoli is defined as 

T\l, m, n) = \l, m, Im + n). (14) 

and all computations are done modulo d. 

Though tangential to the main purpose of this pa- 
per, the above gate-set is another answer to the ques- 
tion posed in Ref. 8|. That is, given a Toffoli, what 
extra gates are required to complete a universal set? Of 
course, the answer provided by the above gate-set in- 
volves measurements in an integral way,and is therefore 
different from the one proposed in Ref. Q . However, the 
above gate-set also addresses the question: Given clas- 
sical computation (i.e. Toffoli and measurements of Z), 
what gates are needed to complete the universal set? 

We now turn our attention to the proof of universality 
for the gate-set presented above. We note that Gottes- 
man has already proven in Ref. that for d prime, ap- 
plying and measuring products of iTs and X J s, plus a 
Toffoli, is universal for quantum computation. All we 
need to do in order to prove universality, is to show that 
we can apply and measure operators of the form X a Z b 
using the above gates. 

Measurements of X followed by measurements of Z can 
produce \i) ancillas for any i. Similarly, we can obtain 
|i) ancillas from measurements of Z followed by mea- 
surements of X. A controlled-sum can be made out of 
a Toffoli by fixing an input to a |1) ancilla. Because a 
controlled-sum is really a controlled- A, fixing the other 
input to |1) produces the X gate. On the other hand, a 
controlled-sum from a state to a |1) ancilla, produces a 
Z on the state: 









X 









The general case of applying X a Z b can be done by a 
series of X and Z gates. All that remains is to construct a 
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method for measuring operators of the form X a Z b . First, 
we note that 

(X a Z b Y = LU abd ( d - 1 )/ 2 x ad Z bd 

{1 d odd , . 

-l ab d=2. (15j 

A. d odd case 

The case d = 2 is rather complicated and will be han- 
dled separately. The general case d odd (remember we 
required d prime) is easy because the eigenvalues of X a Z b 
are the d th roots of unity just like those of X and Z. As 
discussed in the review of phase estimation, being able to 
apply a controlled-A a Z b , combined with measurements 
in the X basis (which includes preparation of X eigen- 
states) is sufficient to measure in the X a Z b basis. 

All that remains is to construct the controllcd-A a Z 6 . 
That is, we need to be able to apply the gate 

|n)®|V>) — |n>® (X a Z b ) n \^) 

= \n) ® x an Z bn u) abn< ~ n - 1 ^ 2 \il>), (16) 

composed of a phase \n,m) — > w '" lm + ab ™(™- 1 )/ 2 | nj TO ) fol- 
lowed by controllcd-sums. The controlled-sum is just a 
Toffoli with an input fixed to one. As for the phase, be- 
cause we have a Toffoli, we have universal classical com- 
putation. We can thus compute brim + abn(n — l)/2 in 
an ancilla, apply a Z to this ancilla, and then erase the 
computation. 

B. d — 2 case 

The d = 2 case is somewhat trickier because our gate- 
set is invariant under complex conjugation, and thus 
there is no way of distinguishing the two eigenstates of 
ZX = iY. We will solve this problem by creating an an- 
cilla that is an eigenstate of ZX, defining it to be the +i 
eigenstate, and then using it to measure and build more 
eigenstates. 

Assume we were given a state 

|*) = ^(|0)+o;|l)) ) (17) 

where u 2 = —1. Clearly, the state is equal to one of the 
two ZX eigenstates: \± Y ) = (|0) ± /\/2. 

Using a controlled-ZA, which is built by the method 
described in the d odd case, we can produce copies of the 
state The idea, similar to the one used for phase 

estimation, is to apply the controlled-ZA from a state 
|0) to the state The target state is an eigenvector of 
ZX with eigenvalue u>, and therefore the relative phase 
is copied over to the first state: 

|6>®|tt) — -^(|0>+o;|l»®|*> 

= |*) ® |*). (18) 



Notice that copying works independently of whether \^>) 
is the +i or —i eigenstate of ZX. Naturally, by subse- 
quently applying a Z, we can also produce the orthogonal 
state |$) = (|0) /V2. 

With our ancilla, we can also measure in this basis. 
This is done by applying a controlled-ZA to the ancilla 
from the state we want to measure: 

I*)® I*) — ► |i}®|*>, (i9) 

I*)® I*) — > |6)<8>|tf), (20) 

and then measuring in the X basis. 

As long as we are consistent in always using the same 
ancilla \*ff), we will have broken the conjugation symme- 
try, and found a way to label, create and measure eigen- 
states of ZX. Of course, we should keep many copies 
of the ancilla, which can be prepared from the original 
state. The operations above also allow us to error correct 
our set of ancillas by copying each, comparing the copies, 
and using majority voting to discard the damaged ancil- 
las. Thus, even if there are some errors in preparation, 
or some of the ancillas decay over time, computation will 
still be feasible. 

All that remains to be explained is how to create the 
first copy of |*). Because a state with a density- matrix 
proportional to the identity can be written as 

P=\l =\\+y){+y\ + \\-y){-y\, (21) 

it is equivalent to having prepared an eigenstate of ZX = 
iY chosen at random. The state p — 1/2 can be produced 
by discarding one qubit of a bell state, and a bell state 
can be produced with a controlled-sum from a |0) an- 
cilla to a |0) ancilla. Therefore, we have shown that we 
can produce the initial eigenstate of ZX, and we have 
completed the proof that the gate-set presented at the 
beginning of this section is universal for quantum com- 
putation. 

IV. UNIVERSAL COMPUTATION FOR SIMPLE 
PERFECT GROUPS 

In this section we will prove that a set of anyons based 
on certain groups can perform universal quantum com- 
putations. Instead of dealing first with the general case 
of non-solvable groups, we will deal with the smaller set 
of groups that are both simple and perfect. 

We remind the reader that non-solvable groups are 
those that contain a perfect subgroup; and a perfect 
group is any non-trivial group, whose commutator sub- 
group equals the full group: [G, G] = G. The property 
of simplicity means that the group has exactly two sub- 
groups that are invariant under conjugation: the triv- 
ial group and the whole group. Because the commuta- 
tor subgroup is invariant under conjugation, it should be 
clear that any simple non-abelian group is perfect. How- 
ever, we shall refer to these groups as simple and perfect 
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to remind the reader that we are dealing with a subcase 
of the general non-solvable case. 

The set of simple perfect groups, which includes the 
groups A n for n > 4, is powerful for computing because 
in some sense we can get from one non-trivial element 
to any other using operations that fix the identity. The 
general case of non-solvable groups will be deferred to 
Section IVII where we will show that a simple perfect 
group can be extracted from a non-solvable group. 

A. Requirements for the physical system 

Here we list the operations, ancillas, and measurements 
that we assume are available on any realistic anyonic sys- 
tem, and which we will use to build our quantum gate-set: 

1. We can braid or exchange any two particles. 

2. We can fuse a pair of anyons and detect whether 
there is a particle left behind or whether they had 
vacuum quantum numbers. 

3. We can produce a pair of anyons in a state that is 
chosen at random from the two particle subspace 
that has vacuum quantum numbers. 

4. We have ancilla pairs of the form \g) ® Iff" 1 ) for 
any g 6 G, where the individual anyons have trivial 
electric charge. 

We remind the reader again that even though all our 
anyons are used in pairs of trivial total flux, we will gen- 
erally only mention one of the anyons of the pair. These 
conventions also apply to ancillas, which means that we 
will refer to the \g) ® Iff" 1 ) state as an ancilla of flux g. 

While the first three requirements are natural opera- 
tions for a laboratory system, it is not clear where the 
ancillas would come from. Depending on the physical re- 
alization there may be many ways of obtaining the ancilla 
reservoir. We discuss one such scheme in Appendix lAl 

B. Computational basis 

Let G be a simple and perfect finite group. Let a and 
b be two non-commuting elements of G. Let d be the 
smallest integer such that a d ba~ d = b. We can assume 
that d is prime, otherwise we could replace a by a d ^ p 
where p is some prime that divides d. 

It turns out that every simple non-abelian group has 
even order. This was first conjectured by Burnside 
|1C| in 1911, and proven by Feit and Thompson |Tlj in 
1963 (in fact, the complete classification of simple finite 
groups was completed in the early 1980's, see for instance 
Ref. 0| )• Using Sylow's theorems, the fact that every 
simple group has even order means that they all include 
a non-trivial element a such that a 2 = 1. Therefore, we 
could always work with a basis of qubits. However, we 
will present the general qudit case both for its elegance, 



and because in some instances a basis of qudits is more 
convenient. 

We will work with a basis of qudits of trivial net flux 

\n) = \a n ba- n ) $ |o n 6- 1 a -n ) (22) 

for < n < d, where we have explicitly described both 
anyons of the pair. 

It should be clear that we can initialize the computer 
by filling up the computational space with |0) ancillas. 
We turn now to the task of constructing the gates pre- 
sented in Section ITTT1 



C. Conjugation by a function 

We begin by describing the technique of conjugation by 
a function, which is especially powerful for simple perfect 
groups. In Scction lll Al we showed that we could perform 
the transformation 

\h)®\g) — > \h)® \hgh~ 1 ), (23) 

where we conjugate the second anyon by the flux of the 
first, while the first anyon remains invariant. We can also 
conjugate an anyon by a product hih-i ■ ■ ■ h n 

Iff) — ► \hth2 ■ ■ ■ Kgh- 1 ■ ■ ■ h^K 1 ), (24) 

where the {hi} are fluxes of other anyons which remain 
unchanged throughout this process. The procedure is 
done by first conjugating by h n , then by h n —i, and pro- 
ceeding leftward until we finally conjugate by hi. 

The above procedure is not terribly useful if all the 
{hi} are fluxes of fixed ancillas, because we could have 
equivalently conjugated by a single ancilla of flux h — 
h\h2---h n . However, some of the fluxes in the prod- 
uct could correspond to anyons that represent qudits of 
unknown state. In this case we can think of the above 
operation as conjugation by a function of the fluxes of 
certain qudits. 

Let us consider what kind of functions can be applied 
in this way. Clearly we are speaking about functions 
that can be written as products of elements of G. The 
elements can include known constants if we use our ancil- 
las to conjugate. We can also include the flux of a qudit, 
which will be of the form a % ba~ % if the qudit is in the 
computational basis (though this may not be the case 
when we are trying to correct leakage). Finally, we can 
include in the product the inverse of the flux of a qudit, 
as discussed in Section Hi Al 

In conclusion, given n qudits with fluxes g\ through g n , 
and a function f(gi, . . . , ffn-i) of the first n — 1 qudits, 
we can conjugate the last qudit by / 

\g n ) — * ■ • - ,9n-i)gnf(gi, ■ ■ -j^n-i) -1 ), (25) 

provided that the function / can be written in product 
form. By product form, we mean that / is a product 
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of the inputs {gi}, their inverses {g" 1 }, and fixed ele- 
ments of G, each of which may appear more than once, 
or not at all. For example, a valid function would be 
/(3i,32) = a 9i^9\ X C 9\ X ^ with a,b,c,d G G. Further- 
more, this transformation does not change the flux of 
the first n — 1 qudits, though it may entangle them with 
the last qudit. 



D. Toffoli gate 

To build the Toffoli gate we must be able to conjugate 
the third qudit by the function / (31,(72), which depends 
on the fluxes of the first two qudits as 



f(a l ba-\a j ba- ] ) = a 1 - 



(26) 



and is arbitrary for values of 31 and 32 that are not in 
the computational basis. If the third qudit is in the state 
k , conjugation by / produces the transformation 



a k ba 



\a k ba- k ) — ► \a ij+k ba- ij - k ), 



(27) 



which is the desired Toffoli gate. 

Given the discussion in the previous subsection, we are 
left with the task of expressing the function / in product 
form. However, it turns out that for simple and perfect 
groups every function has such an expression: 

Theorem: If G is a simple and perfect finite group, 
then any function f(gi,...,g n ) : G n — > G can be ex- 
pressed as a product of the inputs {gi}, their inverses 
{g^ 1 } and fixed elements of G, any of which may appear 
multiple times in the product. 

Not only does the above theorem prove that Toffoli 
gates are possible for any simple and perfect group, but 
it directly proves that any classical function can be com- 
puted. 

The proof of the theorem, which is mostly construc- 
tive, is somewhat long and will be deferred to Appendix 
151 However, to make this seem plausible to the casual 
reader, we would like to illustrate the basic steps needed 
to build a Toffoli gate for qubits. 

The main idea behind the construction is that the func- 
tion / is basically a logical AND of the inputs. A com- 
mutator makes a good logical and because it equals the 
identity if either of its inputs are the identity. Further- 
more, the commutator function can be expanded as a 
product of its inputs. Therefore, we would like the first 
input to take values 1 or c and the second input to take 
values 1 or d, with the requirement that d not commute 
with c, so that we can put them into a commutator. 

Let c/i denote the flux of the first qubit, and 32 the flux 
of the second qubit. Each takes values gi € {b, aba^ 1 }. 
Define the new variables g' l = gib^ 1 E {l,c}, where c = 
[a, 6] = aba~ 1 b~ 1 . It is sufficient to show that we can 
express the Toffoli function as a product of g[, g' 2 , their 
inverses and fixed ancillas. 

Choose an element d that does not commute with c and 
define e = [c,d]. Imagine we could find two functions of 



one element, that can be expressed in product form, such 
that 



hi(c)=d /ii(l) = l, 
/12(e) = a h 2 (l) = l. 



(28) 
(29) 



Using these functions, the Toffoli function can be writ- 
ten as 



/(5i,52) = h 2 



9i, 



(30) 



which when expanded out is a product of the correct 
form. 

The existence of the functions hi, which is discussed 
in more detail in the full proof of the theorem, is a con- 
sequence of G being simple. For any element c S G, the 
group generated by its conjugacy class C(c) is a normal 
subgroup. Because G is simple, this subgroup must equal 
the full group. Therefore, every element d <E G has an 
expression of the form d — x\cx~^ x x 2 cx 2 x ■ ■ ■ XnCX^ 1 for 
some n and some elements {xi} £ G. We can use the 
expression to construct h\: 



Ms) 



x 1 gx 1 1 x 2 gx 2 1 ■ 



' x ng x n 



(31) 



and a similar construction builds h 2 . 

For a concrete example we can work with G = Ac,. 
We begin by choosing an element a, which must satisfy 
a 2 = I, if we wish to work with qubits (d — 2). Because 
of the symmetries of the group, all choices are equivalent 
to a = (12)(34). The next step would involve choosing an 
element b that does not commute with a, and an element 
d that does not commute with c = [a, b]. While any 
choice can produce a Toffoli, the required h\ function 
will be simplified if we can make c and d fall in the same 
conjugacy class. The same can be said for hi if e = [c, d] 
and a are in the same conjugacy class. 

At this point, a little trial and error yield b — (345) 
and d = (234). The computational basis is now defined 
as 



|0) = |6) = |(345)), 

|1) = late- 1 ) = |(435)), 



(32) 



and the remaining group elements are fixed as 

c = (06a- 1 ) = (435)(435) = (345), 

e = (ede" 1 ) d- 1 = (245)(324) = (25)(34). (33) 

The hi functions, which are the only non-constructive 
part of the proof, can be built as simple conjugations 
because of the choices we made earlier: 



hi(g) = h 2 (g) = (521). 9 (125), 



(34) 



where both happen to be the same function by coinci- 
dence. Putting all the steps together we have a function 



/ (3i,52) = (521 



31(435), (521)3a(435)(125) 



(125) 
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(521)ffi(435)(521)ff 2 (435)(125) 
(345) 3 r 1 (521)(345)g 2 - 1 (125)(125) 
(521) 5l (14352) 52 (124) 3 r 1 (15342).g 2 - 1 (521) 
which can be applied with nine elementary conjugations. 



Eventually, one will disappear into the vacuum, and the 
inverse of the flux of that ancilla is the result. Even in 
the presence of errors, this measurement will have a good 
fidelity because the probability of failure is exponentially 
small in the number of fusions. 

A final note is that, because we are always dealing with 
pairs of fluxes, what fusion really means is that we fuse 
the first anyon of our qubit with the first anyon of the 
ancilla. 



E. Measuring Z 

The basic idea behind measuring in the computational 
basis is that if we fuse a flux with another flux of the in- 
verse group element, there is a finite chance that they will 
have vacuum quantum numbers and disappear. On the 
other hand, if the product of the two fluxes is not unity 
then there must be a particle left behind to carry the 
remaining flux (i.e., the total flux is always conserved). 

At this point it might be useful to remind the reader 
why a fusion of g with g~ x will not always turn into the 
vacuum. The short story is that the combined state is 
not invariant when another flux encircles them, implying 
that they have an electric charge component. The state 
that has vacuum quantum numbers is invariant under the 
effect of all fluxes, and hence is the sum of all the states 
in the conjugacy class of g, with the same phase. We 
can figure out the probability of fusion into the vacuum 
by calculating the overlap of the vacuum state with the 
state of two anyons to be fused. The result is 



P=|(Vac(C)|(| ff )®| ff - 1 »| 



1 



\C{g)\ 



(35) 



where C{g) is the conjugacy class of <?, and the vacuum 
state was defined in Section Iff Bl 

Because one fusion will only probabilistically tell us the 
desired result, we should repeat the measurement many 
times to obtain a sufficient degree of accuracy. Besides, 
if we are working with qudits with d > 2 we need to test 
fusion with at least two different fluxes. We therefore 
need to have many copies of the state to be measured. 

Because of the no cloning theorem, copying cannot be 
done exactly, but the transformation 



(36) 



means that we can measure each of the separate copies 
in the Z basis and expect to get the same answer. The 
above transformation can be done with a controlled-sum 
(Toffoli with one input fixed to 1 1)) from the original state 
to a |0) ancilla. Repeating this controlled-sum with many 
target ancillas will produce the above entangled state. 

To summarize, the procedure for measuring in the 
Z basis is first to create an entangled state using a 
controlled-sum. Then try to fuse each of the qudits with 
one of the inverses of the fluxes that are Z eigenstates. 



F. Constructing the zero eigenvector of X 

For the next gates, we are going to need a supply of 
states that are eigenvectors of X with zero eigenvalue: 



d-l 



-y 



(37) 



We will produce them out of pairs of anyons with vac- 
uum quantum numbers. As usual we will just discuss 
one member of the pair, and assume that the equivalent 
operations are being performed on the other anyon. 

One of the possible states that (when paired) have vac- 
uum quantum numbers is the sum of fluxes in the con- 
jugacy class of b. This is approximately what we want. 
Sadly, in general, a state created from the vacuum will 
be a mix of this desired state plus other states, includ- 
ing states that involve dyonic particles (particles with 
both electric and magnetic charge). We will have to fil- 
ter through all this noise to get our X eigenstate. 

The procedure that we will describe below is effectively 
an incomplete swap, that has been extended to the full 
Hilbert space in a logical way. In the computation basis, 
the operations act as 



I*) 



x- 



X 



I*) 



which performs a swap provided that the second qudit 
started in the |0) state. Outside of the computational 
basis, though, the operations are chosen so that we can 
detect whether we obtained the desired |0) state or not. 

We start with two qudit states, one created from the 
vacuum and one which is a |0) ancilla: 



|Vac) <g> |0) = (C|0) + £>|tf j.» <8> |0), 



(38) 



where \^x) is a state orthogonal to the computational 
subspace. If the vacuum pair decohered into a superse- 
lection sector other than the one that contains the com- 
putational basis, the constant C will be zero. This will 
not be a problem as we will be able to detect this case, 
and then start again from this step. 



9 



Using the theorem from Section llV Dl we can conjugate 
the |0) ancilla by a function of the flux of the vacuum pair 
that has the following form: 

f{a l ba- 1 ) = a\ 
/(anything else) = I, (39) 

which is essentially a controlled-sum that has been prop- 
erly defined outside the computational basis. 

The state of the combined system after conjugation 
will be 

. , d — 1 d—1 

-=V \a i ba- i )®\a i ba- i )+Y D t \t> l± }®\a l ba~ % ), (40) 

where {Di} are some constants, and are states 

perpendicular to the computational basis. Note that the 
states \ i Si±) for i ^ are the ones that have flux a % ba~ % 
but have non-trivial charge. The state |^o_l) includes all 
the other fluxes and charges. Depending on the super- 
selection sector in which the vacuum state was created, 
many of the constants C and {D{\ will be zero. 

Now we conjugate by / _1 from the ancilla to the vac- 
uum state yielding 

C\b) ®|6) ®Wba- 1 ), (41) 

i=0 

where {|*q ± >} = {|*o_l)} and the states {\^' i± ),i > 0} 
have flux b but non-trivial charge. 

Now we try to fuse the first qudit with an ancilla of flux 
6 _1 and trivial charge. The only state that can fuse into 
the vacuum with the ancilla is \b), and this will happen 
with finite probability. Note that the ancilla can never 
vanish into the vacuum with a state with charge because 
there is no way of extending the basis to be invariant 
under the stabilizer group of the flux. 

In the end, if the particles disappear into the vacuum, 
the ancilla is left in the desired X eigenstate. Other- 
wise, we repeat the procedure from the beginning until 
eventually the state appears. 

G. Choosing a d th root of unity 

Before we continue building our gate-set, we have to 
address a problem that appears for d > 2, similar to the 
problem that occurred for d = 2 when proving that the 
gate-set is universal. 

So far, we have defined everything in terms of to, a 
non-trivial d th root of unity. But there are d—1 of these, 
and there is a symmetry which interchanges them. We 
will have to break this symmetry by using an ancilla. 

In particular, we need an ancilla that is an eigenstate 
of X with eigenvalue not equal to 1. We will then define 
this state to be the |1) state in the X basis, i.e., 

1 d-1 

I^^E^I*)' ( 42 ) 



which has eigenvalue u, thus fixing our root of unity. We 
then define the other X eigenstates by 

d-1 

|n) = ^=5>-™K>, (43) 

and the operator Z by \n) — ► \n — 1). 

How do we produce the first 1 1) in terms of which ev- 
erything is defined? We start with a |0) (which is al- 
ways well defined and which we know how to construct 
from the previous section), and we apply a controlled- 
X^ 1 (which is a classical function, and thus computable 
from the Toffoli) from this ancilla to a |0) ancilla, which 
produces the output 

|6)®|0)— > -7=£|i>®|i). (44) 

If we discard the second state, we will have a mixed 
state that is a combination of the different X eigenstates. 
This is equivalent to being handed an arbitrarily chosen 
X eigenstate, which we will call 

We can obtain copies of this state by applying a 
controlled-X^ 1 from a |0) ancilla to this state, which 
applies the transformation 

|6) \i) — ► \i) <8> \i). (45) 

We can thus build arbitrarily many copies of the state. 
We still have to worry that this might be the |0) state. 
However, below in the section for measuring X, we will 
give a procedure to detect the |0) which does not rely on 
having |1) ancillas. If we determine that i — 0, we throw 
away all the copies and start over (this will only happen 
with probability 1/d). Otherwise, we relabel our state as 
|1), fixing a value for u>. 

Because we can copy the |1) state, and below we will 
also show how to measure it, we can build a reservoir 
of ancillas in this state, which will be used for all future 
computations. We can even use copying, comparing, and 
majority voting to error correct our reservoir, thus allow- 
ing for computation even in the presence of noise. 

H. Measuring X 

The last gate needed for universality is the measure- 
ment of X. The basic idea is to fuse the pair of anyons 
that form the state to be measured. The |0) eigenstate 
will have some overlap with the vacuum, and will vanish 
with probability p = d/ \C(b)\, where C(b) is the conju- 
gacy class of b. 

The other X eigenstates have zero probability of van- 
ishing because \i) = ^= J2i ^" l \a l ba^ 1 ) is orthogonal to 

the vacuum for i > 0. To detect the state \i) we first ap- 
ply a Z l and then use the above fusion procedure. The Z 
gate can be applied as a controlled-sum with a |1) target 
as discussed in Section ITTT1 
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Of course, the above will require us to have many 
copies on which to measure, which means we need to 
perform the transformation 

J^Cil?) —►5^C i |i}<8|»)®|?}®-- •»!?}, (46) 

i i 

which is done using a controlled-AT -1 with a |0) ancilla 
as control and the state to be copied as target. 

To perform the measurement non-destructively, we can 
fuse all but one of the copies of the state. Alternatively, 
using the Z gate and |0) ancillas, we can always produce 
the rest of the |z) states. The rest of the logic is similar 
to the Z measurement procedure. 

Having completed the construction of the universal 
gates, we have proven that universal quantum compu- 
tation can be performed with anyons from simple and 
perfect finite groups. We now turn to the question of 
whether these operations can be performed in a fault- 
tolerant fashion. 



excitations goes as e~ A / T , where A is the mass gap and 
T is the temperature. 

A good anyonic quantum computer should therefore 
have L » m and T « A. In some implementations, 
however, it may be more practical to accept a small er- 
ror rate from the hardware, and then correct it using 
standard quantum error correction techniques. For such 
cases, we present below the necessary steps needed to 
implement software based error correction for anyons. 

While any of the error correcting codes can be used, 
most techniques require embedding a code space inside 
a Hilbert space on which we can do universal quantum 
computation. However, in the case at hand, our com- 
putational states are embedded in a Hilbert space (the 
states with arbitrary flux and charge) in which we can- 
not perform universal quantum computation. Therefore, 
before starting the recovery protocol, we must first deal 
with states that have leaked out of the computational 
subspace (the subspace in which we can perform univer- 
sal computations). 



V. LEAKAGE CORRECTION 

In this section we will discuss both the motivation and 
the techniques needed to implement error correction and 
fault tolerance in the software of an anyonic computer. 
The main result will be the construction of a leakage 
correction circuit for anyons, which enables the use of 
the standard techniques for handling errors. 

A. Motivation 

Any quantum system that uses non-locality to protect 
its data will be susceptible to errors if a large number of 
its local components are damaged simultaneously. The 
probability for failure is generally exponentially small in 
the size of the system, and is zero in the theoretical limit 
of an infinite system. However, all physical systems are 
finite. Furthermore, practical considerations may force a 
given setup to have a size such that the error of proba- 
bility is small but non-negligible. 

In the case of anyons, errors can occur due to quantum 
tunneling, which is an effect of the high-energy degrees of 
freedom that were frozen out to obtain a two-dimensional 
discrete gauge theory. The probability of this type of 
error goes as e~ mL , where m is the mass of the lightest 
particle that can mediate a charge interaction, and L is 
the separation between anyons. 

Finite temperature effects are another source of error. 
These effects involve the creation of charge pairs from the 
vacuum. Because these pairs have trivial total charge, 
even if they braid with a computational anyon, the net 
charges of the collective three particle excitation will still 
be correct. However, if one of these particles separates 
from the group, or separately braids with another anyon, 
then errors will be introduced. The density of the thermal 



B. Implementation 

To deal with leakage errors we can construct a ver- 
sion of the swap-if-leaked gate described by Kempe et al. 
[T^|. The idea behind the gate is to implement a projec- 
tive measurement that can distinguish the computational 
subspace from its complement. If a state if found to be in 
the computational subspace, it is left alone. Otherwise, it 
is replaced with an arbitrary ancilla that is in the com- 
putational subspace. The ancilla will still be an error, 
but one that is correctable by standard quantum error 
correcting codes. In fact, the general methods of quan- 
tum error correction and fault-tolerant computation can 
be applied to anyons as long as we can reliably project 
leaked qudits into a state in the computational subspace. 

We again focus on the case of simple and perfect 
groups, and defer the general leakage correction protocol 
to the next section. In the current formalism, the compu- 
tational basis is the set of states of a pair of anyons with 
zero total magnetic charge, where each anyon has zero 
electric charge and a magnetic flux of the form a l baT l or 
its inverse. 

The first type of error that we will deal with, is 
when the total magnetic flux of the pair is non-trivial. 
This is a particularly grievous error because, if we drag 
around a pair with a non-trivial net flux, we could be 
introducing errors into all the other qudits. Further- 
more, our assumption that we can perform the operation 
h, g — > h,hgh~ l relied on the fact that the second pair 
had zero net magnetic flux, so it is important that we 
detect and fix this error first. 

To detect a net flux, we take an ancilla \g) ® \g~ x } and 
encircle it by the qudit we are performing the leakage 
correction on. The ancilla will get conjugated by the net 
flux of the qudit, and the qudit will get conjugated by 
the net flux of the ancilla which should be zero. We then 
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fuse the ancilla with a pair with opposite 1 flux. If the net 
flux of the qudit is in the stabilizer of <?, the fusion will 
have vacuum quantum numbers with a finite probability, 
whereas if the conjugation changed the flux of the ancilla, 
there will always be a particle left behind. If we repeat 
this many times with many different ancillas \g) £g> Iff -1 ), 
with good statistical confidence we will be able to tell if 
the net flux of the qudit is in the stabilizer of g. Because 
G has no center, the intersection of all stabilizers is the 
identity, and hence repeating the above with sufficiently 
many different elements g, we can detect a non-zero net 
magnetic charge. 

If we detected a net flux, we replace the state with 
an ancilla in the state |0). Of course, we must be very 
careful when moving the damaged ancilla pair out of the 
region of qudits, so as not to damage other states. That 
is, when moving past other anyons, we always do so in 
the direction in which the damaged pair gets conjugated 
and the good qudits are unaffected. 

In the case when the qudit passes the above test, then 
we have projected into the zero net flux subspacc, but 
otherwise left the state unchanged. The next step is to 
deal with electric charge. Because it is very difficult to 
measure the electric charge of a single anyon, we will start 
with a fresh ancilla |0), made from two anyons neither of 
which have electric charge, and copy the state over. Once 
again we will be using the incomplete swap circuit 



I*) 



10} 



x- 



X 



10} 



I*} 



when acting on the computational basis. Of course, the 
heart of a leakage detection algorithm is how to extend 
the operations outside of the computational subspace. 
The procedure cannot be described simply by a circuit, 
and therefore we will present a way of completing the 
controlled-sum gate so that the above operation will al- 
ways yield a state that is in the computational subspace. 

The following procedure is almost identical to the one 
used to produce |0) states. This is because |0) states 
are obtained by taking a vacuum state and projecting to 
the computational basis, which is primarily leakage de- 
tection. The main difference is that when doing leakage 
detection, we only get one chance of using the qudit (be- 
cause of the no-cloning theorem), but if the state leaked, 
it is acceptable to replace it by anything in the compu- 
tational basis. The latter is clearly not acceptable when 
creating |0) ancillas. 

We will use the incomplete swap procedure for the sec- 
ond round of leakage detection. Recall that by this point 
we have projected the qudit into the zero net flux sub- 
space. Take the qudit and a |0) ancilla, and conjugate 



the ancilla by a function of the qudit's flux: 

f(a l ba- 1 ) = a 1 
/(anything else) = I. 



(47) 



This is the same extension of a controlled sum that was 
used to produce |0) ancillas. 

Afterward, we conjugate the original qudit by /(</) _1 , 
where g is the flux of the ancilla. Note that because 
we know at this point that the original qudit has net 
flux zero, the state of the ancilla will not exit the com- 
putational basis during this operation (though it might 
change within the computational basis if the original 
state had non-zero electric charge). The result of the 
past two controllcd-sums is: 



W®|0> + |V±> 



d-l 

— >|0>®|V||>+X>±i>®|i>, (48) 

8=0 

where parallel and perpendicular refer to inside and out- 
side the computational basis, and none of the ?/> states 
are normalized. Finally, we replace the original pair with 
the ancilla pair and discard the original pair. 

Clearly, the new state will be in the computational 
basis. Furthermore, if the original state was in the com- 
putational basis, then the new state will be equal to the 
old state, and unentangled with the old anyons. 

Having complemented our gate-set with a leakage cor- 
rection scheme, we have proven not only that we can do 
universal quantum computation with anyons, but that 
these computations can be made fault-tolerant. 



VI. UNIVERSAL COMPUTATION FOR 
NON-SOLVABLE GROUPS 

We will now generalize the results of the previous sec- 
tion to any non-solvable group. Unfortunately, in our 
proofs for the simple perfect case, we made extensive use 
of the fact that we can compute any classical function 
simply by multiplying the inputs with ancillas. This is 
no longer true, even if we restrict ourselves just to per- 
fect groups that are not simple. The quickest example is 
^5 x ^5 which is perfect, but has two normal subgroups 
given by each of the A§ factors. Thus, if our two inputs 
are lxl and g x 1, there is no expression made out of 
products in which the results differ in the second factor. 

The above example can easily be fixed by working 
within one subgroup. In general, though, even this is 
not possible, as not all perfect groups have a perfect and 
simple subgroup. However, the following theorem comes 
to the rescue: 

Theorem: If G is a non-solvable finite group, then there 
exists a normal subgroup P of G and a subgroup N, 
normal in P, such that P/N is perfect and simple. 
Once again we defer the proof to the appendix. 
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What the theorem tells us is that we want to work 
with cosets of N in P. That is, we would like to replace 
our old flux eigenstates with states that are labeled by 
elements in P/N and invariant under N. A good guess 
would be 



(49) 




neN 



where x is an element of P/N, and x' is an element in 
the coset that x represents. More specifically, if / : P — > 
P/N is the canonical epimorphism that maps elements 
to cosets, then wc require that f(x') — x. The particular 
choice of x' has no effect on the above definition. 

The above is a good guess but not quite right. A given 
coset may intersect many different conjugacy classes of 
G, each of which lies in a different superselection sector. 
Thus, we are effectively working with mixed states. 

Remembering that we really want to keep our anyons 
in pairs of zero net flux, the right choice for the new 
states is 



1 



Px 



\N\ 



ceC(G) 




(50) 



where again x is an element of P/N, and C{G) is the set 
of conjugacy classes of G. 

These states have the nice property that when conju- 
gated by any element h! S P (or equivalently, when a flux 
h! G P is dragged around them), the effect only depends 
on the coset f(h') of h' , and generates the transformation 



Pf(h')gf(h')- 



(51) 



Because of this, if we use the usual scheme of passing 
one pair of anyons in between another, and they are both 
prepared in states of the above form, the net effect is that 
the inner pair will get conjugated by the outer pair as 



Ph® Pg ► Ph ® Phgh- 1 : 

keeping the pair unentangled. 



(52) 



A. New requirements for the physical system 



distinguished by conjugation, since cgx{cg)~ l = gxg -1 . 
Thus, it may be a difficult problem to distill these flux 
eigenstates from the vacuum. 

A more reasonable assumption is to require the ex- 
istence of ancillas only for the fluxes in the perfect sub- 
group. Another improvement might be to assume that we 
only have ancillas in the mixed states p x defined above, 
where x G P/N. These states might be easier to pro- 
duce because they are obtained from the vacuum by first 
throwing away the anyons with flux not in P or with non- 
trivial charge, and then projecting to a definite coset of 
N in P. Therefore, we will replace our old requirement 
for the existence of flux ancillas by: 

4' We have ancillas in the state p x for any x G P/N . 

It would be highly desirable to be able to prove that re- 
quirements 1 through 4 are sufficient to create the states 
in 4'. Unfortunately, it appears that requirements 1-3 
combined with 4' may neither be a subset nor a superset 
of requirements 1-4. Thus, in a sense, we are imposing a 
different set of requirements for this section. One amelio- 
rating fact is that in the case when P is simple, the states 
p x are just flux eigenstates. We therefore could have used 
requirement 4' for all sections of this paper. We will not 
attempt to describe in the appendix a protocol by which 
these generalized ancillas can be created, however. 



B. Universal computation 



As in Section llVBI we choose two elements a, b G P/N 
such that a d ba~ d = b for some prime d, and a l ba~ l ^ b 
for < i < d. We then define our computational basis 
states as 



Pi — Paiba- 



(53) 



which we define as eigenstates of the Z operator. The X 
operator is defined by the action X{pi) = X piX^ — pi+\ 
and its eigenstates can be obtained using the projection 



operators Pi — Y^j=o u l X i /d by 



d x P l p P, 



(54) 



While the operations of braiding, fusion, and vacuum 
pair creation described in Section liV Al all seem like rea- 
sonable requirements to demand from the physical sys- 
tem, the requirement of flux ancillas is somewhat harder 
to justify. 

In particular, take the case of a group that has a non- 
trivial center, which can occur even if the group is perfect. 
Consider two fluxes g and eg that differ by multiplication 
of an element c in the center. These two fluxes cannot be 



At this point proving that universal quantum compu- 
tation can be achieved is fairly straightforward, and is 
almost identical to the discussion in Section IIVI The 
major differences occur when we have to deal with states 
outside of the computational basis, that is, when creating 
Pq states and when dealing with leakage correction. Both 
of these issues will be dealt with in the next subsection. 
As for the rest of the operations, we will only give a very 
brief discussion: 
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Because the p x states have the same braiding proper- 
ties as those of the fluxes of a group P/N (and in particu- 
lar two Z eigenstates remain unentangled after braiding) , 
the same method for producing a Toffoli applies to them. 

Measuring Z is easy because the pi states have support 
in orthogonal subspaces. The copy (using the Toffoli) 
and fuse with ancillas procedure will work just as well as 
before. 

For the interested reader, we will carry out below some 
of the calculations needed to deal with X eigenstates 
and prove universality. Most of the results seem almost 
miraculous when expressed in the language of density op- 
erators. However, the reader should bear in mind that 
we are only using density operators to account for the 
different superselection sectors. If we just fixed a supers- 
election sector for each particle, we would be dealing with 
pure states, and all of the proofs from the past section 
would carry through. 

We begin by studying the action of the controlled-^ -1 
on X eigenstates: 

-. d-l d-l 

Prn®Pn = ^ E E ^""^"^Vo^ ® Pn 
i=0 j=0 
1 d-l d-l 

- * d E E ^~ im+im ^P^ ® X-'paX-tt 

i=0 j=0 
d-l d-l 

E E u-" n+jm X l Po X^ ® u>- in pnu> jn 

i=0 j=0 
d-l d-l 

»=0 j=0 



which is equivalent to its action on pure states. There- 
fore, once we have p states, we can use the same trick 
as before to break the symmetry and obtain a d th root of 
unity. That is, we do a controlled-^ -1 from the pg with 
a po target to create the state 

Pb ® Pa 

( d-l d-l \ (d-l d-l \ 

i=0 j=Q J \n=0 m=0 / 
1 d-l d-l d-l d-l 

^ 2 E E E E xi p° x " ® u- in+jm p n p pi 

i=0 j=0 n=0m=0 
d-l d-l 

= rf E E (PnPoPl) ® (PnPoPL) (56) 
n— m— 

and then discard (trace out) the first state to get the state 
p = J2 n PnPoPn = J2 n Pn/ d, which gives us an unknown 
eigenstate of X as before. We then discard and repeat 
if we obtained the pg state, and otherwise we relabel the 
state as p\. 



Once the pj state is available, we can use a controlled- 
sum to produce the Z gate, which will allow us to produce 
any X ancilla including more pj states. 

Finally, measuring X works by fusing the pair of 
anyons, because the pi are orthogonal to the vacuum for 
i > 0. The full measurement proceeds as before by copy- 
ing, permuting states using the Z gate, and then fusing. 



C. Leakage detection and generation 

One final issue remains: How do we measure whether a 
state is in the computational subspace? Projecting onto 
the computational subspace is useful because a pg is just 
the projection of a vacuum state to the computational 
basis. Furthermore, this projection will allow us to per- 
form leakage correction. 

One of the new issues that arises for general non- 
solvable groups is that if we have a state in the com- 
putational basis, and we braid it with an electric charge 
carrying a non-trivial representation of the subgroup TV, 
then the state will move outside the computational basis. 
The other issue is that the conjugacy class of an element 
in G might be larger than the conjugacy class of the el- 
ement in P, though given that P is normal, the first set 
will be entirely contained in P. 

Let us begin by examining how the leakage correction 
algorithm must be changed. The first step is to detect 
whether the net flux or charge of the pair of anyons we 
are working on has a non-trivial effect on the states p x . 
The procedure is to braid the pair around the ancilla pair 
and then fuse the ancilla with another ancilla in the state 
p x -i . If the anyon pair has an effect on the ancilla states 
p x , then the fusion statistics will be altered, and this will 
be detectable after many repetitions. If our state is found 
defective we discard it as usual, and replace it by a state 
in the computational basis. Otherwise, we move on to 
the next step. Note that if the anyon pair had a net flux 
in the subgroup N, or in some element outside of P that 
commutes with P, then the state will still advance to the 
next round of error correction. However, this anomalous 
flux or charge will not affect the usual braiding proper- 
ties. 

The second round of error correction is a swap with 
an ancilla in the po state. Note that using our universal 
classical computation in P we can guarantee that if the 
original state was in P, the ancilla ends up in the compu- 
tational basis. However, if the original anyons are outside 
of P, we will get a state that is within P (because P is 
normal) but not necessarily in the computational sub- 
space. The final step is to perform a swap with a second 
ancilla in the po state, where now we know that the first 
ancilla had to be composed of anyons with no charge, 
and fluxes only in P. This guarantees that the final state 
of the second ancilla is in the computational basis, and 
equals the original state if it did not leak, completing the 
leakage correction procedure. 

To create pg we also use a swap, this time between a 
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pair created from the vacuum and a po ancilla. We then 
try to fuse the leftover vacuum state with a p b -i . If they 
fuse into the vacuum, then the ancilla is in a po state. 
The logic is as follows: if the vacuum pair had electric 
charge when created, then the swap will not change the 
charge, and hence it cannot disappear into the vacuum. 
If the vacuum pair has no electric charge but is outside 
of P, then the ancilla is still guaranteed to be in P. Fur- 
thermore, when conjugating the vacuum state, we will be 
conjugating by an element in P. The vacuum state will 
end in a flux state outside of P, which is orthogonal to 
p b -i . Finally, if the vacuum pair is a pair of fluxes in P, 
then it will be of the form pg, possibly superposed with 
other states p x outside the computational basis. But the 
generalized swap can guarantee that a state in P out- 
side of the computational basis, will remain outside of 
the computational basis (just like in the simple perfect 
case). Only when the ancilla is in the state Pq can the 
fusion into the vacuum occur. 

The above procedure for producing pg ancillas com- 
pletes the gate-set for non-solvable groups, and proves 
the main result of this paper: that anyons with fluxes 
in a non-solvable group can perform universal quantum 
computation. 



VII. CONCLUSIONS AND OUTLOOK 

While we have shown that universal quantum compu- 
tation is theoretically feasible for any non-solvable group, 
it is still not yet clear whether we will ever be able to build 
an anyon based computer. First of all there is the fact 
the smallest non-solvable group is A5 which has 60 ele- 
ments. Obtaining such a group from symmetry breaking 
seems problematic. 

One may wonder whether we can do computation with 
solvable groups. For abelian groups, each superselection 
sector consists of just one state, so it is not possible to 
encode quantum data in a topologically invariant fashion. 
Attempts such as Ref. Q encode the quantum data in 
a superposition of position eigenstates, but this has no 
more robustness than using superpositions of positions of 
any other neutral particle of the same mass. 

Hope still remains for solvable non-abelian groups. 
While producing Toffoli gates using conjugation as in this 
paper will most likely no longer be feasible, Toffoli gates 
might still be performed by employing magic states. In 
fact, Kitaev has such a procedure for the group S3 [Tlj . 
The full set of groups which can perform universal quan- 
tum computation remains unknown, but we believe it 
does not include every non-abelian group. 

Furthermore, there are anyons that are not based on 
the electric and magnetic charge model (quantum double 
of a group) presented here. Some of the more exotic 
anyons are likely to be good quantum computers, but in 
general, their computational power remains unknown. 

We have also neglected to present in this paper an 
account of the resources used to perform computations. 



While it should be clear that computations can be done 
with at worst a polynomial overhead in the size of the 
input, some gates (in particular those that require calcu- 
lations of arbitrary functions over the group) may require 
resources that are exponential in the size of the group. A 
lot of the wasted resources may come from the description 
in terms of general groups, though. For a fixed group, the 
resources can probably be significantly reduced. 

Finally, there remains the question of physical systems 
which contain anyons. Because of the requirement of two 
dimensions, we must look for quasi-particles in some two 
dimensional medium. There are some indications that 
non-abelian anyons may arise in the fractional quantum 
Hall effect (see Refs. 0,0] and references therein). How- 
ever, at the moment, there are no physical systems out 
of which the anyonic computer may be built. Even if 
no physical implementations are ever found, though, this 
subject will hopefully still be interesting because of its 
beautiful mix of computation, particle physics, and group 
theory. 



Acknowledgments 

The idea for universal classical computation with sim- 
ple and perfect groups was initially suggested by Alexei 
Kitaev, to whom I am highly grateful. The author would 
also like to thank John Preskill, Jim Harrington, and 
James Chakan for proofreading this paper, and Meg 
Wessling both for translating this paper into proper En- 
glish and for many useless but entertaining conversations. 

This work was supported in part by the National Sci- 
ence Foundation under grant number EIA-0086038 and 
by the Department of Energy under grant number DE- 
FG03-92-ER40701. 



APPENDIX A: CREATING THE ANCILLAS 

As discussed in the main text, the requirement of a 
supply of calibrated flux ancillas needs further justifica- 
tion. In this section we will show that for a perfect and 
simple group, the requirements of braiding, fusion, and 
vacuum pair creation can be supplemented by one extra 
measurement to allow the distillation of flux ancillas. We 
will not cover the general non-solvable case, though. 

The new measurement involves determining whether a 
single anyon has trivial flux or not. Indeed, this measure- 
ment may even be done destructively. The plausibility of 
this measurement relies on the fact that non-zero flux 
charges are topologically non-trivial configurations that 
often have much higher masses than their electric charge 
counterparts. Naturally, dyons also have large masses 
and will be detected as having non-trivial flux. 
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Step 1: Creating electric ancilla pairs 

The procedure for creating flux ancillas begins by cre- 
ating single anyons with zero flux. These are obtained 
by creating a vacuum pair, measuring the flux of the 
first particle of the pair, and discarding the second one if 
the first one had non-trivial flux. 

The next step is to create pairs of anyons, where each 
anyon has zero flux and unknown charge, but the pair 
has vacuum quantum numbers. Of course, if we could 
non-destructively distinguish trivial from nontrivial flux, 
we could skip this step, as the vacuum pairs always have 
vacuum quantum numbers. 

Take two of the single electric charges we have pro- 
duced. We are going to try to project this state onto the 
desired state with vacuum quantum numbers. Consider 
the process of creating a pair of anyons from the vac- 
uum, braiding one of them around the pair of charges, 
and then fusing the vacuum pair. If the pair of charges 
had vacuum quantum numbers, then the vacuum pair 
will remain in the vacuum state throughout this process, 
and fuse into the vacuum at the end with unit probabil- 
ity. On the other hand, if the pair of charges does not 
have vacuum quantum numbers, then there will be a fi- 
nite probability that the pair created from the vacuum 
will leave a particle behind after fusion (since the vac- 
uum is the only state that is left invariant by the action 
of every flux) . 

Repeated application of this process will be a projec- 
tive measurement which determines whether the pair of 
charges has vacuum quantum numbers. If we project 
onto a vacuum pair, then we have found a good charge 
ancilla pair. If the pair does not project onto the vac- 
uum state (because the two anyons do not transform in 
conjugate representations, or because we projected to a 
state orthogonal to the vacuum), then we pair them up 
with other charges and repeat the process. While slow, 
this process will eventually yield as many electric charge 
pairs with vacuum quantum numbers as needed. 



Step 2: Identifying the magnetic charges 

The electric ancilla pairs are useful because they can 
perform a non-destructive measurement of magnetic flux. 
The procedure is to take a member of the electric charge 
pair, drag it around the anyons or group of anyons whose 
total flux we want to measure, and then fuse it with its 
pair. 

To describe the effect of the fluxes, let R(g) be the rep- 
resentation of the first electric charge of the pair. Let |n) 
be an orthonormal basis for the space on which R acts, 
and let be the dual basis for the conjugate repre- 
sentation R* under which the second charge transforms. 
The effect of a flux g is then 

£>)®|n*)^X;(iJG7)|n>)®|n*). (Al) 




FIG. 3: Using electric charges to check if gi = gi. 



Just as before, if the total flux is non-trivial, there will 
be a good chance that the fusion of the electric charges 
will leave a particle behind. On the other hand, if the 
total flux is trivial, even if the total charge is not, the 
pair of electric charges will always fuse into the vacuum. 

Repeated application of this procedure will determine 
whether the total flux is trivial or not. Furthermore, this 
procedure will at worst introduce decoherence in the flux 
basis, but will leave all flux eigenstates unchanged. 

We can use this procedure to compare the fluxes of 
two anyons. In particular, consider two pairs created 
from the vacuum. Measure the total flux of the first 
anyon of the first pair combined with the second anyon 
of the second pair. If the combined flux is trivial, the 
first anyon of each pair has the same flux; otherwise the 
flux is different. The procedure is depicted in Figure 

The above procedure allows us to sort the flux pairs 
into "bins" that depend on the total flux of the first anyon 
of the pair. We will get as many bins as elements of 
G, each containing an unlimited supply of vacuum pairs 
which carry the same flux in the first anyon of the pair. 
At this point, if the fluxes have not decohered in the flux 
basis, then we must have an entangled state involving all 
anyons in a given bin. Throwing away a single flux from 
each bin will produce the desired decoherence, just as it 
did when breaking the various symmetries in the main 
part of this paper. 

All that remains is to identify each bin with an element 
of G. Assume that we were given an assignment of an ele- 
ment of G to each bin. The assignment could be checked 
by using the following procedure. First, we note that 
any finite group G may be described by a set of elements 
{gi} and a set of relations of the form gi t ■ ■ ■ gi n = 1 which 
they obey. To check that the assignment is correct, we 
just need to check all the relations (supplemented by the 
trivial one element relations gfg™ — g™ m )- These can be 
checked again with the electric charge ancillas, using a 
loop that circles each of the fluxes in the relation in the 
correct sequence. 

To generate guesses, we could just randomly assign to 
each bin an element of g, which gives us a probability of 
success of at least 1/(|G|)L Of course, we can be a lot 
smarter, as the above procedure can help us figure out the 
powers of a given element (including its inverse) and even 
the elements in its conjugacy class. Thus the need for 
guesswork is minimal, and some of the choices correspond 
to different valid assignments (i.e., automorphisms) of the 
group. 
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Analysis of the produced ancillas 

At this point we have almost produced the desired an- 
cillas, with one caveat: the individual anyons do not have 
trivial charge (i.e., they may be dyons). However, all we 
have done to the pairs, after creating them from the vac- 
uum, is to measure the flux of one of the anyons. That 
means that the electric charge portion of the state is still 
in the vacuum state. More technically, if the ancilla pair 
circles a flux that commutes with the flux of one anyon 
of the ancilla, then the state remains unaltered. This is 
the same behavior that the pure magnetic charges would 
have. 

Some careful thought at this point shows that these 
states are good enough for the quantum computation 
procedure presented in the bulk of the paper. Indeed, 
going back and repeating all the steps with these gener- 
alized ancillas would require very few modifications. The 
fusion to measure in the Z basis would now have a lower 
success probability, which is compensated by a higher 
rate of producing acceptable |0) states, but otherwise 
most gates remain unaltered. We have therefore suc- 
ceeded in constructing an ancilla reservoir, which, while 
slightly different then the one initially desired, is useful 
for universal quantum computation. 



APPENDIX B: MATHEMATICALIA 

This appendix proves the major mathematical theo- 
rems needed in the bulk of the paper. We begin by stat- 
ing the definitions of some of the mathematical terms 
used: 

Perfect group: A non-trivial group G such that [G, G] = 
G. Note that [G, G] is not the set of elements of the form 
[51,32] = 5i525r 1 <?2~ 1 but rather the group generated by 
these elements. Even if G is perfect, there may not be a 
commutator expression for every element. 
Non-solvable group: A group that has a perfect subgroup. 
Normal subgroup: A subgroup H of a group G such that 
ghg~ x £ H for every h £ H and g £ G. 
Simple group: A group with no normal subgroups other 
than the whole group and the trivial group. 

Before we get to our main theorem, we will prove a the- 
orem that will allow us to deal with general non-solvable 
groups. We intend to show that we can extract from non- 
solvable groups a simple and perfect group. The simple 
perfect groups (which can also be described as the simple 
non-abelian groups) are the ones on which we can per- 
form universal classical computation and are therefore 
important for this paper. 

We begin by defining the n th derived subgroups by the 
relations = [G^-V, G^'V] and G^ = [G,G]. A 
solvable group is one for which G w = {1} for some 1. 
A non-solvable group must have an i such that for every 
j > i, GW = G^> and C?W is non-trivial. The group 
C?W is perfect, thus the definition for solvable groups 



is consistent with the definition for non-solvable groups 
given above. 

Furthermore, all the groups 

G (n) 

are normal subgroups 
of G. This can be proven by recalling the property 
g[gi,92}g~ 1 = [ggi9~ 1 ,gg2g~ 1 ]- The rest follows by in- 
duction because G^ is normal in G, and G^' is normal 
in G if G^ -1 ^ is. We have therefore shown that every 
non-solvable group G has a perfect normal subgroup P. 

Sadly, this subgroup is not necessarily simple. How- 
ever, we can prove that every perfect group P has a nor- 
mal subgroup N such that P/N is perfect and simple. 
We choose N to be a normal proper subgroup of P such 
that no other normal proper subgroup of P has more el- 
ements, which is well defined because P is finite. Let 
/ be the canonical epimorphism P — > P/N which maps 
elements into cosets. Because / is surjective we have 
[P/N, P/N] = [f(P),f(P)\ = f([P,P}) = f(P) = P/N, 
which, combined with the fact that P/N is non-trivial, 
shows that P/N is perfect. 

Finally, assume that P/N has a normal, non-trivial 
proper subgroup A. Then B = f^ 1 (A) is a normal 
subgroup of P, because for any elements b\ , 62 £ B 
and p £ P, we have /(61&2) = fibijffa) £ A and 
fipbiP' 1 ) = f{p)f(bi)f(p)- 1 £ A. Furthermore, B is 
a proper subgroup of P, and N — / _1 (1) is smaller than 
B = f^ 1 (A), leading to a contradiction. Therefore, P/N 
is simple, and we have finished proving the following the- 
orem: 

Theorem: If G is a non-solvable finite group, then there 
exists a normal subgroup P of G and a subgroup N , 
normal in P, such that P/N is perfect and simple. 
Proof: Shown by the above text. 

We now turn our attention to using our groups to com- 
pute classical functions. We shall prove that the set of 
functions that can be written in product form is com- 
plete, in the sense that it includes every function from 
G™ — > G, if G is simple and perfect (or equivalently 
simple and non-abelian). This was first proven in the 
mathematical literature by Maurer in 1965 ^(|. In the 
computer science literature, a related result was proven 
by Barrington |17|. In this paper, we will provide our 
own constructive proof for the following theorem: 
Theorem: If G is a simple and perfect finite group, then 
any function f(gi, . ■ . , 3,1) : G™ — > G can be expressed as 
a product of the inputs {g{\, their inverses {g" 1 } and 
fixed elements of G, any of which may appear multiple 
times in the product. 

Proof: Throughout this proof we will refer to the set 
of functions that can be expressed in the above form as 
"computable." Proving the above statement is equiva- 
lent to showing that all functions are computable. The 
proof consists of building a series of computable delta 
functions that map most elements to the identity, and 
then expressing arbitrary functions as a product of these 
delta functions. 

Step 1: Given a group element a not equal to the iden- 
tity, let G(a) denote its conjugacy class. Then the sub- 
group generated by the elements of C(a) is equal to G. 
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This is because the subgroup is a nontrivial, normal sub- 
group of G and G is simple. 

Step 2: Fix two disjoint subsets A and B of G. Define a 
family of functions {8^' B {g) : AU B — > G} with elements 
labeled by c € G: 



?A,B 



(a) 



5?' B (b) 



VaeA 
VbeB. 



(Bl) 



If the function S^' B is computable for some 1, then 
every function in the family is computable. To prove this 
choose any d 6 G. By Step 1 there is an expression for d 
as a product of elements in the conjugacy class of c (for 
instance, d — gxcg^gicg^c). Then 5 d ' B is obtained by 
substituting 5f' B for c in the expression. 
Step 3: Fix a set A, an element b not in A, and an 
element x ^ b. If a function <5^ ,B is computable for 
some B such that b G £?, then there exists a computable 
function ^ ,B with two new sets such that A\J {x} C A' 
and b E B' . The function can be obtained from 



(.9) = [<*d' S (ff))^ 



(B2) 



using Step 2. The above equation assumes that we have 

A B 

extended the domain of S d ' to G, which can be done in a 
natural way once we have fixed a product representation 

A B 

for 6 d ' . The element d was chosen to not commute with 
bx^ 1 . Such an element must exist because G is simple 
and non-abelian, and hence has no center. The element 
e is just e = [d, bx^ 1 ]. 

Step 4: The functions defined by 8 b c (g) = 5^' B , with 
A = G — {b} and B = {b}, are computable. To prove 
this start with A\ = {1} and B\ = {b}. The function 



§M,Bi j s computable because it is in the same family as 
/(flO — .9 — Sg' Bl - Then proceed by induction, using 
Step 3, on the elements in G — {b} that are not included 
in Ai. 

Step 5: For a fixed set of ordered elements b\, . . . ,bi 
define a family of functions labeled by c: 



,gi) = c g! = bi,..., and g t = b t 
,gi) = \ otherwise. (B3) 



The same proof in Step 2 shows that if any function of the 
family with c 7^ 1 is computable, then the entire family 
is computable. 

Step 6: Fix i € Z + and elements 61,... € G. If 

the function #c (<?i, ■ ■ ■ ,gi) is computable, then so is 



^...fc 

sufficient to be able to compute 



the function 5c 1 '"°* +1 (ffi) • • • iSitfl'i+i)- By Step 5 it is 



rbi...i>i. 



^ 1 - 6 *(fli,...,ft),C +1 (fli+i)J . 

(B4) 

where the function S d ' +1 (gi+i) is computable by Step 4, 

and d is chosen so that e = [c, d] 7^ 1. 

Step 7: Using induction on the number of inputs of 

the function, and starting from the base case S bl (gi), 

it is clear that all the functions defined in Step 5 are 

computable. 

Step 8: Every function is computable because 



/ (pi, • • • , = n " ■ n 6 )Xbu...M) s<) • 

Q.E.D. 
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